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Abstract 



JiT I The objective of this paper is to discuss the dynamical instabiUty 

Q^ ■ in the context of Newtonian and post Newtonian regimes. For this 

t^^ . purpose, we consider non-viscous heat conducting charged isotropic 

fluid as a collapsing matter with cylindrical symmetry. Darmois junc- 
^Ij ' tion conditions are formulated. The perturbation scheme is applied 

^\| . to investigate the influence of dissipation and electromagnetic field on 

the dynamical instability. We conclude that the adiabatic index F has 

smaller value for such a fluid in cylindrically symmetric than isotropic 

sphere. 
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1 Introduction 

The stability problem of relativistic stars has been an interesting subject for 
researchers. The static stellar model would be interesting if it remains stable 
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under fluctuations. This can be well discussed in terms of adiabatic index F. 
The dynamical stability PP-|1] of isotropic spheres in general is given by 

4 M 
r > - + n— , 
- 3 R' 

where n is a number of order unity that depends on the structure of the star 
and M and R are the mass and radius of the star respectively. For example, 
for white dwarfs n = 2.25. 

Different physical aspects of fluids play a key role in the dynamical in- 
stability and evolution of self-gravitating systems. It was shown [5] that 
the dissipation increases the instability range at Newtonian corrections but 
makes the fluid less unstable at relativistic corrections. Chan et al. [B] found 
that anisotropy and radiation affect the instability range at Newtonian and 
post Newtonian (pN) regimes. The same authors [7J explored that the shear- 
ing viscosity decreases the instability range of the fluid both at Newtonian 
and pN regimes. Chan [8j studied collapsing radiating star with shear viscos- 
ity and concluded that shear viscosity would increase anisotropy of pressure 
as well as effective adiabatic index. Dev and Gleiser [S] investigated the dy- 
namical instability of neutral fluid sphere with varying energy density and 
showed that anisotropy enhances the stability whenever the tangential pres- 
sure is greater than the radial pressure. Sharif and Kausar [10] discussed the 
dynamical instability of expansionfree fluid at Newtonian and pN order in 
f{R) gravity. 

In a strong gravitational field, a star requires more charge to be in stable 
form. The effect of electric charge in a self-gravitating bound system has 
been studied by Rosseland [11]. Bekenstein [12] found that charge density 
has a significant effect on the system but unable to show whether the star is 
stable or not. Glazer [13] investigated the effect of electric charge upon the 
dynamical instability of isotropic fluid. Zhang et al. [H] showed that the 
electric charge changes the structure of neutron star when the charge and 
mass density are same. Further, the presence of electric field increases the 
stability of the star [T5] . 

Esculpi and Aloma [IB] found that both anisotropy and electric charge 
would boost the stability of sphere under radial adiabatic perturbation as 
compared to isotropic distribution of matter. Ernesto and Simeone [ITj in- 
vestigated that the inclusion of charge extends the stability regions for both 
bubbles and shells around a black hole. Joshi and Malafarina [18] studied the 
stability of Oppenheimer, Snyder and Datt black holes under small tangential 



pressure. They found that this black hole is not stable within the collapsing 
cloud whose final fate of collapse is a naked singularity rather than a black 
hole. In a recent paper [19], we have discussed the effects of electromagnetic 
field on the dynamical instability of a spherically symmetric expansionfree 
gravitational collapse. 

Some recent work [20]- [23] indicate great interest in cylindrical gravita- 
tional collapse with different fluids with and without electromagnetic fleld. In 
this paper, we use cylindrically symmetric distribution of collapsing isotropic 
fluid to discuss the dynamical instability with electromagnetic fleld. Darmois 
junction conditions [21] are used to match the interior non-rotating cylindri- 
cally symmetric spacetime to exterior cylindrically symmetric spacetime in 
the retarded time coordinates. The paper is organized as follows. In section 
2, we formulate the Einstein-Maxwell fleld equations. We discuss junction 
conditions in section 3. In section 4, the perturbation scheme is applied on 
the fleld as well on the dynamical equations. Section 5 provides the discus- 
sion of the dynamical instability in terms of non-dissipative and dissipative 
perturbation. We conclude the results in the last section. 

2 Interior Fluid Distribution and the Field 
Equations 

We consider a collapsing cylinder bounded by a hypersurface S. The line 
element for the interior region has the following form 

ds^_ = -A^dt" + B^dr^ + C^dO^ + dz\ (1) 

where A, B, and C are functions of t and r. To preserve the cylindrical 
symmetry, we have the following constraints on the coordinates 

— oo<t<oo, < r < oo, — oo < z < oo, 0<^< 27r, 

where we assume comoving coordinates inside the hypersurface. The interior 
coordinates are taken as {x"} = {^i ^; ^; z}{a = 0, 1, 2, 3). The source in the 
fleld equations is assumed to be locally dissipative in terms of heat flow given 
as [5], [25] 

'^afS = (a* + P)UaUi3 + PQa/S + QaUfi + QfiUa, (2) 

where /i is the energy density, p is the isotropic pressure, «„ and q^ are 
four-velocity and radial heat flux vector of the fluid respectively satisfying 



<iaU°' = 0. Moreover, we have 

u" = A-^d^, g" = q6'^, u'^Ua = -1. (3) 

The Maxwell equations can be written as 

Fa/3 = h,a - 0a,/3, F;;" = 4nJ", (4) 

where (pa is the four potential and J" is the four current. The electromagnetic 
energy-momentum tensor is given by 

Ea^ = ^ (f2F^, - ^F-y'F.sg.p^ , (5) 

where F^^ is the Maxwell field tensor. Since the charge per unit length of 
the cylinder is at rest with respect to comoving coordinates, the magnetic 
field will be zero in this local coordinate system [23], [26]. Thus we can write 

0, = (0(t,r), 0,0,0), J" = C«", (6) 

where (p{t, r) and ( is the scalar potential and charge density. The conserva- 
tion of charge requires 

s(r) = 47r /" CBCdr (7) 

Jo 

which is the electric charge interior to radius r. Using Eq.([T]), the Maxwell 

equations (jl]) yield 

^0' = AnCAB', (8) 

/ = 0. (9) 

Here dot and prime represent derivatives with respect to t and r respectively. 
Integration of Eq.® with respect to r, assuming 0'(t, 0) = 0, gives 

<P' = —s{r). (10) 

The electric field intensity is defined as 

E{t,r) = ^,. (11) 
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The Einstein-Maxwell field equations, G^a = Sn (T^^ + E^^j , for the 
interior metric gives 
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B^C E'^C AB^C A^BC 
Equations (fT5|) and (fT6|l yield 

A2 yXc 'EC) ^ 52 \Ac ~ ~^c) ^ c \B^ ~aJ~' ^ ^ 

Thorne [27] defined C-energy for the cylindrical symmetric spacetime 

m(t, r) = E{t, r) = -(1 - r'^V^fVaf), (18) 

8 

The circumference radius p, specific length / and areal radius f can be defined 



as 



p' = e(l)a6"l) , l^ = ^{2)a^(2) , ^ = pi 



(19) 



where ^(1) = ^, ^(2) = ^ are the Killing vectors and E{t,r) represents the 
gravitational energy per unit specific length of the cylinder. The specific 
energy of the cylinder [28] in the interior region with electromagnetic field 
can be written as follows 



m{t,r) = E'{t,r) 



1 



1 + 



C 



c 

E 



s 
2C' 



(20) 



The dynamical equations can be obtained from the contracted Bianchi iden- 
tities (T~"^ + E~"^).^jj = yielding the following two equations 

f, + q'A+(^2^ + ^ + ^^qA+{^, + p)(^ + ^] = 0,(21) 



3 Junction Conditions 

In this section, we take a timelike 3-D hypersurface E which splits two AD 
manifolds V~ and V^ as interior and exterior regions respectively. The 
interior region is defined in Eq.(Il]) while the exterior region is described by 
cylindrically symmetric manifold |29] in the retarded time coordinate as 

\ R R^ J 

where 7^ = — g-, A being the cosmological constant and M^u) is the total 
mass, Qiy) is the charge inside the boundary surface S and v is the retarded 
time. For smooth matching of the interior and exterior regions, Darmois 
conditions ED can be stated as follows: 



1. The continuity of the line elements over S 

{dsl)^={dsl)^={ds^)^. (24) 

2. The continuity of the extrinsic curvature over S 

[K,,] = KJ-K^=0, (z,j = 0,2,3). (25) 

The boundary surface S in terms of interior and exterior coordinates can be 
defined as 

/-(t,r) = r-rs = 0, (26) 

U{v,R) = R-R{vj,) = Q, (27) 

where r^ is a constant. Using Eqs. 026|) and f l27|) . the interior and exterior 
metrics take the following form over S 

[dsV)^ = -A\t,r^)dt^ + C\t,rj:)de^ + dz\ (28) 



(ds^, 



2M{u) Q^{iy) 



R{i^h R{^\ 



2^^^]du' + RUde' + ^'dz'). 



dv 



The first condition implies 

df 1 

- = A{t,T^)~\ C{t,r^) = R^{v) = -, 
ar 7 

yielding the radius of S in both interior and exterior coordinates 

\2 
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2M{u) ^ Q{vf ^ JR^ 



Ry: ' i?s ' du 
The outward unit normals to S are found by using Eqs. fl26|) and ([2 
n- = (0,S(t,rs),0,0), 

' dR^ 
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(31) 

(32) 
(33) 



The non- vanishing components of extrinsic curvature in terms of interior and 
exterior coordinates are 
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(35) 



dR\ f dv 
dr J \dT 
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(36) 



Using the continuity of extrinsic curvature and Eqs. (l30l) and (I3T!) . we get 

2 
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(37) 



(38) 



where s = Q has been used. This shows that the difference of both masses 
at the boundary E turns out to be | which is zero in spherical case. This is 
the least satisfactory condition of the C-energy where the isotropic pressure 
is equal to the radial heat-fiow. It is mentioned here that pressure at the 

Y] 

boundary will be zero only if there is no dissipation, i.e., q = 0. 



4 The Perturbation Scheme 



A stable stationary black hole solution under perturbation indicates the final 
state of the dynamical evolution of a gravitating system. This scheme is also 
used to identify basic properties of black hole solution [30] ■ Here we perturb 
the field equations, Bianchi identities and all the material functions upto first 
order in A such that the fiuid is in hydrostatic equilibrium form. These are 
given by [5] -[7] 



Ait,r 


) = Ao(r) + AT(t)a(r), 


(39) 


B{t,r 


) = Bo{r) + XT{t)b{r), 


(40) 


Cit,r 


) = rB{t,r)[l + XT{t)c{r)], 


(41) 


E{t,r 


) = Eo{r) + XT{t)e{r), 


(42) 


fi{t,r 


) = l^o{r) + XHt,r), 


(43) 


p{t,r 


) = po{r) + Xp{t,r), 


(44) 


q{t,r 


) = Ag(t,r), 


(45) 


m{t,r 


) = mo(r) + Am(t,r), 


(46) 



where < A ^ 1. The static configuration of Eqs. flT2|) - flT4j) and (H 
obtained by using Eqs. (l5^ - (H51) as follows 
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The corresponding perturbed field equations become 
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Bianchi identities ( ETj) and (12^ for static and perturbed configuration yields 

1 5' 
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pd + (/^o + Po)-r 
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(54) 



An 



-^ + ^"j gAo + ("^ + c ) (/io + Po)r = 0,(55) 
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Using Eg. 0511) in flS^ and then integrating, it follows 

1 Ao 
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Similarly, the static configuration of Eq. (120 p turn out to be 
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The above equation gives the total energy entrapped inside S. From Eqs. flM]) 
and fl59l). we have 
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The second order linear differential equation can be obtained from Eqs.(l5T 
and (152|) with the help of junction condition and the relation po = 



T + 2I3T -aT = 0, 
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The general solution of Eq. (l62p is given by 



T{t) = - exp[(-/3s + V«s + /?|)t], 



(63) 



i;)t^^) 



(65) 



where for the sake of real solution and range of instability, we assume that 
«£ > and /3s < 0. We see that the system starts collapsing at t = — oo with 
T(— oo) = in static position and it continues to collapse with the increase 
oft. 
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5 Dynamical Instability 

Here we discuss the dynamical instability at Newtonian and post-Newtonian 
regimes. We relate F in terms of p and Ji as the ratio of specific heats 



p = r — ; — ^l, 



(66) 



r is taken to be constant throughout the dynamical instability. Using Eq.([5 
in the above equation, we obtain 



P = -Po^ {— + c]T- —-^T ■ ( 



{r)T. 



(67) 



Also, from Eq. (l52p . it follows that 
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Making use of Eqs.( !5T|) . ( 1571) . ( 1671) and ( l68|) in ( !56l) . we get the collapse 
equation 



26 \ 1 Ao^ Po 

-'^^^Vo^V-^.^^H^TV.' 

8n B^ ^0 + Po 
/26 
V-^o y '" " fj-o + Po 



26 
B'o 



+ c) (/io + Po 



|,(„„,,,(i,|)-' 



X 



in(^ + c] TpoB'o - r-^^Ao0(r) + -^ 






V-DQ 



A'of b 



c - Sn^EoBoieBo + bEo) + 167r6SoPo 



Ao \B, 

+ (Eocy+(Eoc + e)(- + ^)+^{rbEoy 

\r BoJ rBo 



Eo[e' 



- i^|4 - ' («^' 



Now we discuss the dynamical instability at Newtonian and pN regimes in 
the following cases of vanishing and non-vanishing heat flow. 
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5.1 Non-Dissipative Perturbation 

In this case, we assume that the fluid is non-dissipative, i.e., q = 0. As a 
result, Eq. flST]) can be written as 
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This is trivially satisfied if both the terms are taken to be identically zero, 
i.e., 
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b = AoBo 
while the remaining two terms yield 



a'-°- 



(71) 



(72) 



c = — 



Pb 



Ei + EqE'o + M\E^ 



\c- csjl - 8mo + WiT^E^Bor. (73) 



,yUo+Po f^o+Po 

Since there is no dissipation, thus Eqs. (!58l) and (!64l) gives 0(r) = and 
/3(r) = respectively. Now using these results and Eqs. (172l) and ( I73ll in the 
collapse equation ( l69ll . we obtain 
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Newtonian limit 



+ Jl-8mo + 167r2S25or 



0. (74) 



For the instability condition at Newtonian limit, we take Aq = 1, Bq = 1, 
/^o ^ Po and neglecting the relativistic effects. Equation (17^ yields 

- 2p[,r + 2p[, - 3(^2 + SttV^o' + ^o^o) + (1 - 87rVE2)as/io = 0, (75) 

Thus the instability condition for the isotropic fluid in the cylindrically sym- 
metric spacetime is given by 



r < 1 + 
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^ I Po I 



(76) 



where we assume p'^ < 0. This equation shows that the electromagnetic fleld 
increases the instability of the collapsing fluid. 



Post-Newtonian limit 

In this limit, we use ^40 = 1 — ^^, Bq = 1 — ^^^ and the relativistic effects 
upto order ^^. Thus we have 

-2p[,r + 2p'o - 3{El + SttV^o" + EoE',) + (1 - S^rE^ [Wnpofio 
X (1 - r) - 87rVo^o(e + ^o) +(1 + c)asfio] - Eq [e' + e(l + Sn^rEl) 
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Since the system starts collapsing as t— ;■ — oo, i.e., T < 0, it follows from 
Eq.( l5T|) that c > 0. Hence the instability condition at pN limit becomes 

r < 1 + -^ [2{El + SvrV^o" + ^o^o) + S^Voe^^o' + ^o (e' 
^ I Po I 

+ e(l + SttV^o') + Kc)] - TT^ (167rpo^o + (1 + c)as/io) • (78) 

^ I Po I 
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This shows that the instabihty range is increased by the first three terms in 
the square brackets and is decreased by the last term. It is mentioned here 
that the first term in the square brackets is the Newtonian term while the 
remaining terms are the relativistic corrections. 

5.2 Dissipative Perturbation 

When g 7^ 0, the perturbed solution corresponding to Eq. fl72|) becomes 

b{r) = AoBo [1 + vfir)] , (79) 

where ?7 > is of the order ^^a. Using the above solution in Eqs.f lSSp and 
flM|) . (f){r) and /?(r) turn out to be of order ^^. 

Newtonian limit 

Applying the same approach as defined in the non-dissipative case and using 
the values Aq = 1, Bq = 1, fio ^ Po, Eq.f l69|) at Newtonian limit reduces to 

-2(1 + r]f)p',T + 2(1 + r//)p[, - 3(1 + r]f)iE', + E^E', + SvrVi^o') 

nf 
+ (l-87rVE2)(l + r//)/ioas + 7^«E = 0. (80) 

In general, the heat flow g > and T < for a collapsing fluid, implying c < 
and /' < from Eq. (l5T]) . Consequently the instability condition becomes 



r < 1 + 
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2b'ol 87r|p^|(l+r7|/|)J 

We note that the electromagnetic field as well as dissipation enhance the 
instability range, which supports the result for spherical case [5]. 

Post-Newtonian limit 

Here we take ^40 = 1 — ^^, i?o = l — ^^ and consider the relativistic 
correction terms upto order ^^a 'jj-^g collapse equation at pN limit yields 

-(2 + 277/ + c)p;r + (2 + 277/ + c)p'^ - (3 + 877/ + c){El + E^E'^ 

+ 87rV^o') + ^^/'«s - (1 - S-K^rEl) [STrVo^oe + 87r(2 + 277/ + c) 

X poUqV - (1 + 77/ + c)as/Uo - 87r(l + 77/)(2po - 7r^o)yUo + {.vf + c') 
X {-p'^ + El + ^0^0 + S^'^^o )] - ^0 [e' + (1 + Sn'^rEl) (e + E^c) 
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Using the same arguments as given in Newtonian limit, we have the instabihty 
condition 
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This imphes that the range of instabihty is increased by Newtonian term and 
relativistic terms in the square brackets while it is diminished by the last four 
terms. 



6 Summary 

This paper has been addressed to discuss the dynamical instability of cylin- 
drically symmetric non-viscous heat conducting isotropic fluid with electro- 
magnetic field. For spherically symmetric spacetime [Ij, it has been shown 
that when a star is compressed or expanded adiabatically, its dynamical in- 
stability depends on the numerical value of the adiabatic index, i.e., |. If 
r > |, the pressure in the star is strong enough than the weight of the lay- 
ers which makes the star stable. On the other hand. If F < |, the weight 
increases very fast than the pressure and the star collapses resulting a dy- 
namical instability. 

We have found that for the cylindrically symmetric case, the dynamical 
instability depends on the critical value 1. When F > 1, the stability of the 
system is obtained while F < 1 yields a dynamical instability. Further, we 
have studied the role played by the electromagnetic field and dissipation in 
the early stages of dynamical instability. It is interesting to mention here 
that different ranges of instability may lead to different patterns of evolution 
of stars. We note from Eqs.l ffGl) and ( ISTl) that electromagnetic field as well 
as dissipation make the system more unstable at Newtonian regime. This 
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supports the fact that different physical aspects of the fluid has a great 
relevance in the evolution as well as structure formation of self-gravitating 
objects. A similar remark applies for Eqs. flTSj) and fl83|) at pN regime. It 
is mentioned here that only relativistic effects are taken into account at pN 
regime. We have found that the results of dissipative case reduce to non- 
dissipative if heat flow vanishes. 
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